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The sectoral synchronization observed for the Japanese business cycle in the Indices of 
CO ' Industrial Production data is an example of synchronization. The stability of this synchro- 

nization under a shock, e.g., fluctuation of supply or demand, is a matter of interest in 
physics and economics. We consider an economic system made up of industry sectors and 
goods markets in order to analyze the sectoral synchronization observed for the Japanese 
' business cycle. A coupled oscillator model that exhibits synchronization is developed based 

(-H , on the Kuramoto model with inertia by adding goods markets, and analytic solutions of the 

stationary state and the coupling strength are obtained. We simulate the effects on syn- 
chronization of a sectoral shock for systems with different price elasticities and the coupling 
^H' strengths. Synchronization is reproduced as an equilibrium solution in a nearest neighbor 

graph. Analysis of the order parameters shows that the synchronization is stable for a fi- 
nite elasticity, whereas the synchronization is broken and the oscillators behave like a giant 
oscillator with a certain frequency additional to the common frequency for zero elasticity. 

> , 

^ ■ §1. Introduction 

\ The business cycle is an example of synchronization, and this has been studied 

in nonlinear physics. In particular, the stability of synchronization under a shock, 
e.g., a fluctuation in supply, is a matter of interest. Sectoral synchronization was 
observed for the Japanese business cycle in the Indices of Industrial Production data 
from 1988 to 2007.M3 By synchronization, we means there is a constant phase 
difference between industry sectors. For instance, if random noise components are 
^ ' removed in an appropriate manner, the business cycle is written as 



Xi = siniiot + 6i), (I'l) 

^ = f, (1-2) 

where Xj, 9i, u, and T are the normalized growth rate of production for sector i, 
the phase of the business cycle for sector i, the common angular frequency, and the 
common period of the business cycle, respectively. 

We consider an economic system made up of industry sectors and goods markets. 
Various types of shock may occur in the economic system of interest. If production 
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in a specific sector suddenly decreases significantly, the imbalance of demand and 
supply may destroy the synchronization. However, if the price of the goods in the 
market is quickly adjusted, a sudden change of production in the sector will be 
absorbed by consumers as a decrease in demand in coupled sectors. 

In this paper, we describe a synchronized coupled oscillator model with flucuat- 
ing goods markets. The effects on synchronization with respect to the magnitude of 
shock, price flexibility, and network topology are studied using the coupled oscillator 
model. The paper is organized as follows. In section 2, existing theories are briefly 
reviewed. In section 3, we explain the formulation of our coupled oscillator model, 
and we describe our analysis of sectoral synchronization in section 4. Finally, section 
5 presents our conclusions. 

§2. Existing Theories 

In this section, existing theories for the synchronization of economic systems, 
power systems, and physical and biological systems are reviewed. It is noted that the 
first two are large-scale manmade systems that exhibit similar dynamical behavior 
in terms of synchronization. 

2.1. Economic Systems 

The business cycle is observed in most of industrialized economies. Economists 
have studied this phenomenon by means of mathematical models, including various 
kinds of linear non-linear, and coupled oscillator models. 

Interdependence, or coupling, between industries in the business cycle has been 
studied for more than half a century. A study of the linkages between markets and 
industries using nonlinear difference equations suggests a dynamical coupling among 
industries.® A nonlinear oscillator model of the business cycle was then developed 
using a nonlinear accelerator as the generation mechanism.'^' In this paper, we stress 
the necessity of nonlinearity because linear models are unable to reproduce sustained 
cyclical behavior, and tend to either die out or diverge to infinity. 

However, it is noted that a simple linear economic model, based on ordinary 
economic principles, optimization behavior, and rational expectations, can produce 
cyclical behavior much like that found in business cycles.^' An important question 
aside from synchronization in the business cycle is whether sectoral or aggregate 
shocks are responsible for the observed cycle. This question has been empirically 
examined, and it was clarified that business cycle fluctuations are caused by small 
sectoral shocks, rather than by large common shocks.EJ 

As the third category of model, coupled oscillators were developed in order to 
study noisy oscillating processes like national economies .'^^ Simulations and em- 
pirical analyses showed that synchronization between the business cycles of different 
countries is consistent with such mode-locking behavior. Along this line of approach, 
a nonlinear mode-locking mechanism was further studied that described a synchro- 
nized business cycle between different industrial sectors.-^ 
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2.2. Power Systems 

Power systems made up of many synchronous machines and transmission hnes 
exhibits synchronization. The stabihty and control of this synchronization of the 
machines has been studied in power system engineering.fiSli 

The dynamical property of a single synchronous machine is often studied in the 
single machine connected to the infinite bus shown in Fig. [I](a), where Vg, Z5, I, xi, 
ri, Voo, and ZO are the generator voltage, generator phase, current, transmission line 
reactance, transmission line resistance, infinite bus voltage, and infinite bus phase, 
respectively. Here the infinite bus is characterized by constant voltage and constant 
frequency. The active power and reactive power Qe are calculated by 

Pe + iQe = Vgexp{i6)I\ (2-1) 

where * indicates the complex conjugate, and the current / is given by 

^ _ Vgexp{i6) - Voo ^2.2) 
n + ixi 

By substituting Eq. (|2-2|) into Eq. (|2-1|) . the active power Pg is obtained as the real 
part of the equation 

p = Ya}^ sin 6 = P"*"^ sin 5, (2-3) 
XI ^ ^ 

by considering relation r/ <C x^. Pg is plotted as a function of the generator phase 5 
in Fig. [1] (b). It should be noted that Pe is the electric power transmitted through 
the transmission line and works as a synchronizing force in power systems. 
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Fig. 1. Single Machine connected to Infinite Bus and Synchronizing Force 



2.3. Physical and Biological Systems 

Many collective synchronization phenomena are known in physical and biological 
systems.!^ Physical examples include clocks hanging on a wall, an array of lasers, 
microwave oscillators, and Josephson junctions. Biological examples include syn- 
chronously flashing fireflies, networks of pacemaker cells in the heart, and metabolic 
synchrony in yeast cell suspensions. 

Kuramoto proposed a coupled oscillator model to explain this rich variety of 
synchronization phenomena.'^^K^''^ In his model, the dynamics of the oscillators 
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are governed by 

N 

Bi = uji + J2 sm(ej - Oi), (2-4) 

where 9i, Ui, and kji are the oscihator phase, the natural frequency, and the couphng 
strength, respectively. The second term of the RHS of Eq. ()2-4p is identical to Eq. 
(j2-3|) If the coupling strength kij exceeds a certain threshold, the system exhibits 
synchronization. 

§3. Coupled Oscillator Model 

In this section, we first explain the formulation of our coupled oscillator model, 
and then give an illustrative example to understand the basic behavior of the model. 

3.1. Formulation 

Our model is developed based on the Kuramoto model with inertia^'^' by adding 
goods markets. We consider the system to consist of oscillator i and oscillator j. 
The angle 6i of oscillator i is written as 

ei = ujt + ei, (3-1) 

where uj and 9i are the frequency and phase, respectively. 

The energy dissipated as heat from oscillator i at a rate proportional to the 
square of the angular velocity is 

Pd = KoiOif, (3-2) 

and the kinetic energy accumulated in oscillator i at a rate proportional to the square 
of the angular velocity is 

p. = (3.3) 

where Ko and / are a dissipation constant and a moment of inertia, respectively. 
From Eq. (j3Tp . the angular difference AOji is written as the phase difference A9ji = 
9j — 6i = 6j — Oi using the phase 6i. The power transmitted from one oscillator to 
another is given by 

Pi = -%sinZ\%, (3-4) 

using the phase difference A9ji. In Eq. ()3-4p . the negative sign is used to indicate 
that power is lost from oscillator i. 

By substituting relations (|3-2p . (j3-3p . and (|3-4p into the power balance equation 
for oscillator i (Pg = Pd + Pa + Pt), we obtain an equation corresponding to Eq. ()2-4p 

Ioj6i = Ps- Kdu? - 2KD0j6i + kji sin AOji, (3-5) 

using the approximate relation 9i <^ to. From the above discussion, without loss of 
generality, we obtain an equation describing the dynamics of an A^-oscillator system, 

N 

h = Pi- a§i + sin (3-6) 
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where Ioj, Ps — Kdu'^, and 2Kj:)L0 are replaced by 1, Pi, and a, respectively. Pi is 
regarded as the net input to oscillator i, which is equal to the difference between 
the input to oscillator i from outside the A^-oscillator system and the output from 
oscillator i to outside the A^-oscillator system. Hereafter 9i is written as 9i for 
simplicity. 

It is noted that the synchronizing force is interpreted as being the goods markets. 
The meaning of the goods markets is clarified by adding the sectoral fluctuations of 
demand or supply, 6ji, as 

N 

Oi = Pi- aOi + ^{kji sin AOji + 5ji}. (3-7) 

i=i 

The introduction of these sectoral fluctuations is intended to be consistent with the 
importance of the small sectoral shocks found in Ref.® Figure [5] depicts a system 
consisting of six oscillators. Figure [2] (a) is the Nearest Neighbor (NN) graph, and 
Fig- E] (b) is the Complete (C) graph. Oscillators and goods markets are indicated 
by circles and rectangles, respectively. In Fig. [2] (a), the number of the markets is 
six, because they are only open to the neighboring oscillators. In Fig. [2](b), however, 
the number of markets is five for each oscillator and the total number of markets is 
equal to (6 x 5)/2 = 15. 



(a) eo(t) (b) Go(t) 




63(1) 63(1) 



Fig. 2. Network Topology 
Demand dij for, and supply Sji of, good i can be written as 

dij = do + 5ij, (3-8) 

Sji = So + (3-9) 

and are determined through the goods market ij. do and sq are the equilibrated 
demand and equilibrated supply of good i in goods market ji, and 6ij and 6ji are 
the sectoral fluctuations of demand and supply, respectively. In Eqs. (|3-8p and p-9p 
the equilibrated demand do and equilibrated supply sq are given by 



do = So = kji\ sm{9* - e*)\, 



(3-10) 
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where the equihbrated phase 9*{i = 1, • • • , A^) are solution of Eq. (|3-7|) with 9i = 
6i = 0. The relations between demand d or supply s and the price p of good i are 
written using the price elasticity of demand or the price elasticity of supply e^, 



d 
do 



(3-11) 



(3-12) 



p_ 

So \Po, 

The price as a function of demand and supply is depicted in Fig. [3l Figure [3] (a) 
is plotted for small value of e^. In this case, the demand is shown as a vertical line. 
Oscillator j responds to the sectoral fluctuations of supply 5ji according to price p. 
Depending on the market flexibility, oscillator j changes its demand by 5ij, given by 



-6ji {ea < 0) 
{ed = 0), 



(3-13) 



by responding to the fluctuation of supply 6ji . 




2 3 4 

Quantity 



Fig. 3. Demand and Supply 



3.2. Illustrative example 

We consider an oscillator system in order to understand the basic behavior of 
the coupled oscillator model. For the NN graph, the model equations are written as 

01 = Pi - aOi + {k2i sm{e2 - Oi) + ,521 }, 

02 = P2- Oi02 + {^21 sin(6'i - 6*2) + ^12} + {A;32 sin(6l3 - 6^2) + 532}, 

(3-14) 

0N^\ = Pn-1 - a6N-i + {kN--iN-2 sm{9N-2 - On-i) + Sn-2N-i}, 
+ {kNN-1 sm{eN - On~i) + 5nn-i}, 
On = Pn - oiOm + {kNN-i sin(6'Ar_i - 6m) + 6nn-i}- 

Analytic solutions of the stationary state will be obtained by solving the simultaneous 
equations of (|3-14p with 0j = 0j = 0(i = 1, • • • ,A^). We consider the case, where 
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Pi = P = 1, Pi = 0(i = 2, ••• ,Ar - 1), Ptv = -P = -1, and On = 0. The 
analytic solutions obtained for the stationary state 9* and the synchronizing coupling 
strengths kji are 

^* = 0*^^ + arcsin-^, (3-15) 




(3-16) 



We simulate the behavior of the system by solving Eq. ()3-14p numerically with 
the initial condition Oi{0) = 9i{0) = 0(i = I,-- - ,N). Synchronization was repro- 
duced as an equilibrium solution in a simple NN graph as shown in Fig. HI The 
response to a sectoral fluctuation at t = 100 is also simulated in the NN graph. In 
the case of zero elasticity (e^ = 0), the synchronization was broken, as shown in Fig. 
[5] (a). In contrast, in the case of finite elasticity (e^ < 0), stability was restored after 
a shift of phase, as shown in Fig. [5] (b). 




50 100 150 200 

Time step 



Fig. 4. Stationary Solution 
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Fig. 5. Fluctuation and Elasticity 

§4. Analysis of sectoral synchronization 

In this section, we analyze the sectoral synchronization observed in the Japanese 
business cycle using the coupled oscillator model described in the previous section. 
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The effects on synchronization resulting from a sectoral shock are studied for systems 
with different price elasticities and coupling strengths. 

4.1. Synchronization 

For the business cycle with a common period of T = 60 months, the phases 
of production, shipment, and inventory are given in the article.'^ The sectors are 
rearranged in decreasing order of the observed phase of production. If we assume 
that the network is the NN graph, we can calibrated the coupling strengths kij using 
Eq. ()3-16p with the observed phases 6* . Analytical solutions of the stationary phases 
are then obtained using Eq. (|3-15|) . 

The analytical solutions of stationary phases are compared with the observed 
phases in Fig. [6] (a). The agreement between the analytical solutions and the 
observed phases is quite good. We simulate the behavior of the system by solving 
Eqs. ()3T4p numerically with the initial condition 9i{0) = ^i(O) = 0(i = 1, • • • ,N). 
Synchronization was reproduced as an equilibrium solution in a simple NN graph, 
as shown in Fig. [6] (b). 




00 0.3 10 1.5 ■ 50 100 150 200 

Observed phases Time step 

Fig. 6. Parameter calibration and synchronization 

4.2. Fluctuating Goods Market 
The complex order parameter 

1 ^ 

q{t) = J^Y1 = r(cos(</.) + isin{cl))) (4-1) 

is defined as a macroscopic quantity that corresponds to the centroid of the phases 
of oscillators.^^' The radius r measures the coherence and 4> is the average phase. If 
Re(g(i)) ~ 1 and Im(g(t)) w 0, the oscillators remain in the synchronization region 
where the phase differences are fairly small. In contrast, if Re(g(t)) and lm(q{t)) 
oscillate between 1 and —1, the oscillators behave like a giant oscillator with a 
frequency additional to the common frequency lo. 

First, we simulate the response to the sectoral fluctuation at t = 100 in the NN 
graph. The result for a finite elasticity (e^; < 0) is shown in Fig. [7] (a). It is seen 
that Ke{q{t)) f« 1 and Im(g(i)) ss 0, even after the sectoral shock was applied in the 
middle of the network at t = 100. This means that the synchronization is stable, 
i.e., the oscillators remain in the region where phase differences are fairly small. For 
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the case of zero elasticity (e^ = 0), Re{q{t)) and lm{q{t)) oscillate rapidly between 
1 and —1 as shown in Fig. [7] (b). This means that synchronization is broken and 
the oscillators behave like a giant oscillator with a high frequency additional to the 
common frequency cj. 



(a) (b) 
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Fig. 7. Price elasticity and stability of synchronization 

Next, we simulate the oscillator system with weaker coupling strengths kij. We 
expect that the system to behave as if the oscillators were uncoupled for kij below 
than a certain threshold kfj. For the Kuramoto oscillator, the exact formula of 
the critical coupling strength kfj has been derived and verified with the results of 
numerical simulations. The results of our simulation with weaker coupling strengths 
of O.Skij are shown in Fig. [HI The time evolution of the phases 6i shown in Fig. [8] 
(a) depicts a few oscillators separating from the main part of the coupled oscillator 
system, which still exhibits synchronization. The oscillator system acts like a giant 
oscillator with a low frequency additional to the common frequency to, as seen in 
Fig. [8](b). The phases 6i at t = 200 in Fig. [8](c) show that the system disintegrated 
into four parts and that the main part moves with small phase differences. 

This result is quite different from our expectation, but this is also reasonable, 
because the system under consideration has Pi = 1, Pj = 0(i = 2, • • • , — 1), P/v = 
— l(iV = 21). Hence, one end of the system is pulled in the positive direction 
and the other end is pulled in the opposite direction. Therefore, if the coupling 
strengths are weak enough, the oscillators in both ends of the system are separated 
from the main part. The oscillators in the main part lose tension and shrink to the 
small phase difference. However, it is noted that the system in the C graph with 
Pi 7^ 0(i = 1, • • • , A^) might behave similarly to the Kuramoto oscillator. 

§5. Conclusions 

The sectoral synchronization observed for the Japanese business cycle in the 
Indices of Industrial Production data is an example of synchronization. The stability 
of this synchronization under a shock, e.g., fluctuation of supply or demand, is a 
matter of interest in physics and economics. 

First, existing theories for synchronization for economic systems, power systems, 
and physical and biological systems were reviewed. 
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Fig. 8. Case with the coupling strength multiphed by 0.3 



We then considered an economic system consisting of industry sectors and goods 
markets. A coupled oscillator model exhibiting synchronization was developed based 
on the Kuramoto model with inertia by adding goods markets, and analytic solutions 
of the stationary state and the coupling strength were obtained. 

Finally, we analyzed the sectoral synchronization observed in the Japanese busi- 
ness cycle using the coupled oscillator model. Effects on synchronization resulting 
from a sectoral shock were simulated for systems with different price elasticities and 
coupling strengths. The agreement between the analytical solutions and the observed 
phases was quite good. Synchronization was reproduced as an equilibrium solution 
in a simple NN graph. Analysis of the order parameters showed that the synchro- 
nization was stable for a finite elasticity (e^ < 0), whereas the synchronization was 
broken and the oscillators behaved like a giant oscillator with a frequency additional 
to the common frequency for zero elasticity (e^ = 0). 

In the future work, we intend to study 

1. Simulation of the oscillator system for the complete graph with Pj ^ 0(i = 

2. Theory of phase transition for synchronization, including derivation of the exact 
formula for the critical coupling strength k^-. 

3. Development of the economic implications for sectoral fluctuations and the 
propagation of risk through the economic network. 
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